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Binary Mixtures of Particles with Different Diffusivities Demix 
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The influence of size differences, shape, mass and persistent motion on phase separation in binary 
mixtures has been intensively studied. Here we focus on the exclusive role of diffusivity differences 
in binary mixtures of equal-sized particles. We find an effective attraction between the less diffusive 
particles, which are essentially caged in the surrounding species with the higher diffusion constant. 
This effect leads to phase separation for systems above a critical size: A single close-packed cluster 
made up of the less diffusive species emerges. Experiments for testing of our predictions are outlined. 

PACS numbers: 61.43.Hv, 64.75.Xc, 64.75.Yz, 82.70.Dd 


Binary mixtures consisting of non-adhesive particles in 
thermal equilibrium can demix if their constituents dif¬ 
fer in size m- This demixing is explained in terms of 
an increase in the entropy of the system. However, non¬ 
adhesive binary mixtures that are far from thermal equi¬ 
librium may demix even when the particles are of equal 
size. Examples include mixtures of active and Brownian 
particles pun], and shaken granular media consisting of 
particles of different masses [mill]. While the impact of 
activity on the phase separation dynamics has been fre¬ 
quently addressed PIlo], the specific role of differences 
in the diffusion constants has not yet been investigated. 

Here, we study conditions under which diffusivity dif¬ 
ferences among non-adhesive spherical particles of equal 
size and shape can lead to demixing. We investigate bi¬ 
nary mixtures of two particle species which differ only 
with respect to their diffusion constants, and interact by 
a short-ranged repulsive force. We find that demixing 
is promoted for large relative differences in the diffusion 
constants and high overall packing fractions. The binary 
mixture exhibits phase separation into a solid-like clus¬ 
ter of the species with the lower diffusion constant (‘cold’ 
particles) and a gaseous phase of the ‘hot’ particles with 
the higher diffusivity; see Fig. As discussed below, 
these predictions may be tested using various experimen¬ 
tal systems including mixtures of granular disks Hffln] 
or mixtures of photo-activated colloids dsj m]. 

The Brownian dynamics of the binary system of A^hot 
hot and A^coid cold particles is described by a set of cou¬ 
pled Langevin equations: dtVi = 

where the sum runs over all particles i S A^} 

with N = A^hot + A/oid- All particles have the same 
mobility ^ but differ in the value of their diffusion con¬ 
stants, Di G {T^hot,T^coid}, which determine the respec¬ 
tive magnitude of the spatially isotropic and Gaussian 
white noise {'nia{t)0jp{t')) = 2Di5ij5ap6{t — t'). The 
interactions between the spherical particles of radius 
a are taken as short-ranged harmonic repulsive forces, 
Fij = k{2a — rij)vij if particles overlap < 2a), 
and Fij — 0 otherwise. Here, k denotes the spring 
constant and = jr^ — r^l the inter-particle distance 
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FIG. 1. Snapshots of particle configurations in a binary mix¬ 
ture of 500 hot particles (orange, light gray) and 500 cold 
particles (blue, dark gray) with diffusion constants that differ 
by a factor D = I5coid/Tlhot = 10“^ for a packing fraction 
4 > = 0.2 at (a) early times t = 6 • lO'^a^/Dhot, and (b) late 
times t = 6.4 • 10“^ a? jD^^ot of the coarsening dynamics. Ini¬ 
tially, small shape-less clusters of cold particles form, which 
later coalesce and merge into a single cluster. Inset: Close- 
up of cluster with color coded local hexatic order parame¬ 
ter 4'6,i = lA/ir^^ with Mi denoting the Voronoi 

neighbors of particle i and 6ij is the ‘bond’-angle between 
particles i and j. 

{iij = (r^ — Yj)/rij). In order to mimic hard particles, 
we use large values of fc • /x, thus ensuring that particle 
overlaps decay quickly (see Supplemental Material [H])- 
Recently, a similar model has been used to study clus¬ 
tering near a hard wall |18] . For specificity, we consider 
A^hot = A/om; we have checked that changing the ratio 
of hot to cold particles by a factor of two in either di¬ 
rection does not lead to qualitatively different results. 
As control parameters we consider the ratio of diffu¬ 
sion constants D := I?coid/T^hot and the packing frac¬ 
tion (j) •= where L is the system size of the 

simulation box with periodic boundary conditions. 

Our Brownian dynamics simulations show that for suf¬ 
ficiently large relative differences in the diffusion con¬ 
stants an initially random distribution of hot and cold 
particles spontaneously segregates into a solid-like phase 
of cold particles surrounded by a gas-like phase of mainly 
hot particles; see Fig. [^a,b) and Supplemental Material 
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for a video m This clustering instability is consistent 
with analytic predictions based on a low density expan¬ 
sion for the same system [19] , In our simulations the 
phase separation process begins with the formation of 
small, shapeless clusters of cold particles of varying size. 
We observe coarsening dynamics as these clusters both 
coalesce and continue to incorporate individual particles, 
finally leading to a single, spherical cluster at large time 
scales. The bulk of this cluster exhibits hexagonal or¬ 
der [inset of Fig. [^b) and Supplemental Material m 
Fig. S8]. Moreover, we have measured the number of cold 
particles in the largest cluster. Moo, at time scales where 
the system has reached a steady state. Figurej^shows the 
fraction of cold particles in the largest cluster at steady 
state, Moo/iVcoidj as a function of D for four packing frac¬ 
tions (j). We find that for small D almost all cold particles 
end up in the largest cluster. However, with increasing 
D we find a pronounced drop in Moo/iVcoid, which sets in 
at a point that strongly depends on the packing fraction 
(j). This suggests that there is a threshold value IIcrit(0)j 
which marks a phase transition from an isotropic phase, 
where cold and hot particles are homogeneously mixed, 
to a demixed phase, in which they are segregated into 
a solid phase of cold and a gaseous phase of mainly hot 
particles. 

To further quantify this phase separation scenario we 
explored the dynamics of cluster formation and growth. 
Particle motions are determined by intrinsic diffusion 
(with diffusion constant Hhot or Hcoid) and particle col¬ 
lisions. For I?hot ^ Dcold^ movement of cold particles is 
primarily driven by collisions with hot particles, lead- 
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FIG. 3. (a) Time evolution of the fraction of cold particles 
in the largest cluster M/Ncom in a system of Wot = Woid = 
300 hot and cold particles with diffusion constant ratio D = 
1.5 X 10~^ at total packing fraction (j) = 0.2. Each color 
(gray shade) depicts a single realization of initially randomly 
distributed particles. The lag time ti and the time for cluster 
growth Tg is shown for the rightmost curve. The dashed line 
indicates ten collision times (b) Cluster size distributions 
of cold particles for systems with (j) = 0.2 and D = 0. Different 
symbols and colors correspond to different particle numbers; 
see legend of (c). (c) Time traces of M/Wow (same system 
as (b)). Black solid line: Numerical solution of Eq. Q with 
parameters as in Supplemental Material |17| . Fig. S2, started 
after the corresponding lag time. 



FIG. 2. Fraction of cold particles in the largest cluster of the 
system, averaged after saturation Moo/Woid (symbols) as a 
function of the diffusion constant ratio D = Dcoid/Dhot of 
cold and hot particles for systems with Woid = Wot = 300 
particles. Different colors correspond to different packing 
fractions (j). Clustering is also observed for packing frac¬ 
tions higher as depicted here (see Supplemental Material m 
Fig. S6). Values Mcxo/Woid > 0.8 (shaded region) indicate 
that a single large and stable cluster of cold particles has 
formed. The saturation value M^a was obtained from 10 re¬ 
alizations. Solid lines show solutions of Eq. 0 for t —>■ 00 , 
where the parameter Qint was fitted to a single set of D and 
(j> (indicated by *); see main text and Supplemental Mate¬ 
rial [T7| . 


ing to an effective mesoscopic diffusion constant Dqs 
of the cold particles. Our simulations show that for 
a single cold particle immersed in a bath of hot parti¬ 
cles at packing fraction i^hot, the effective diffusion con¬ 
stant can be approximated by Zleff = ^(t^hotDdot with 
7 « 0.28; see Supplemental Material [T7], Fig. SI. As il¬ 
lustrated in Fig. I^a), four qualitatively different regimes 
can be discerned for the phase separation process. Ini¬ 
tially, on a time scale determined by the collision time 
Tc = L^/(4Woid.C)eff), there is fast assembly of many 
small clusters of cold particles. These small clusters co¬ 
exist for an extended lag time which typically is much 
longer than the collision time Tc. Subsequently, for the 
system sizes considered, we observe coarsening dynam¬ 
ics which finally leads to the formation of a single cluster 
whose growth saturates due to depletion of cold particles. 
While the duration of the growth period Tg is similar for 
different realizations of the stochastic dynamics, the lag 
time Tf, varies significantly. 

The large variance in T(, suggests that the coarsen¬ 
ing dynamics is preceded by a cluster nucleation phase, 
i.e. a cluster must first reach a critical size before it 
can stably grow. As this implies that no stable clus- 
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ters can emerge in finite systems with particle numbers 
close to or below the number required to form a viable 
nucleus, we measured the weighted cluster size distribu¬ 
tion p(m) = m • rim/Ncoid {nm- frequency of clusters 
of size m of cold particles) for different particle num¬ 
bers N while keeping the parameters cj) = 0.2, D — 0 
and iVhot/-^coid = 1 fixed. We hnd that for small sys¬ 
tems (iVcoid = 50), there is a broad distribution of cluster 
sizes originating from the continuous assembly and disas¬ 
sembly of clusters [Fig. ib, c)], while for larger systems 
(iVcoid = 150), a single and stable cluster can form. We 
conclude that coarsening indeed requires the formation of 
a critical nucleus whose size is of the order of A^coid ~ 150 
for the parameters considered in [Fig. ib,c)]. 

After nucleation, growth of the cluster is driven by 
a balance of attachment and detachment processes. To 
quantify these processes we develop a phenomenological 
mean-held theory for the dynamics of the cluster mass 
M{t). Assume that the cluster is approximately spher¬ 
ical with radius R and hexagonally packed with pack¬ 
ing fraction rj = 0.9 [see Fig. Ha), inset] such that 
Ma^ « rj-KR? . In two dimensions, the hux of diffusing 
particles towards a sphere of radius R scales as the diffu¬ 
sion constant times the area density of the particles and 
is independent of R [20]. Hence the rate of accretion of 
cold particles with effective diffusion constant D^g scales 
with Watt(A'/) = £>eff(A^coid - M)/L'^, where A^coid - M 
is the number of cold particles outside the cluster. De¬ 
tachment may be mediated either through collisions of 
hot particles with the cluster surface or through intrinsic 
diffusion of cold particles. Similarly to the case of cold 
particle attachment, the detachment rate due to colli¬ 
sions with hot particles scales with WcoU = DhotNgot /L'^, 
i.e. the flux of hot particles towards the cluster. For 
the detachment rate due to intrinsic diffusion we take 
Wi„t(M) = N„/{(l>hotTint), where Tint = a^/Dcoid is the 
time it takes a cold particle to leave the cluster due to 
its intrinsic motion and Ncr = TrRja « /r] is the 

number of cold particles on the cluster surface. More¬ 
over, because hot particles block cold particles from dif¬ 
fusing away from the cluster we multiply the correspond¬ 
ing time-scale by i^hot to account for this caging effect. 

The three processes combined lead to the time evolu¬ 
tion of the cluster mass: 

dtM = Oiatt ^att (Af) Ocoll ^coll Oint a^int(Af) , (1) 

where aatt, ctcoii and aint denote dimensionless coeffi¬ 
cients, which were determined by fitting the solution 
of Eq. 0 to time-traces M{t) obtained from ten inde¬ 
pendent realizations of our Brownian dynamics simula¬ 
tion for a single parameter set; see Supplemental Ma¬ 
terial HZ], Fig. S2 for details. We obtain ctatt = 82.5, 
acoii = 0.03 and aint = 7.5 • 10“^. Using these values we 
find that the solutions to Eq. 0 are in good agreement 
with the simulation time-traces for different packing frac¬ 
tions or cold particle diffusion constants [Fig. [^c) and 


Fig. S2 in the Supplemental Material mi Overall, this 
analysis shows that cluster growth is limited by collision- 
mediated diffusion of cold particles, and that cold surface 
particles are strongly attached to the cluster structure, 
with detachment events being rare. Moreover, the mean- 
field analysis also gives the correct saturation values for 
the cluster mass Mqo in the regime where a single and 
stable cluster develops [Fig.l^ shaded area]. As expected, 
outside of this regime Eq. Q deviates from the simula¬ 
tions, since it does not take the coexistence of multiple 
clusters of fluctuating size and shape into account. A 
complete theory of cluster growth would need to incor¬ 
porate cluster fragmentation and coalescence as well as 
Ostwald ripening [2T] . 

The observed saturation of cluster growth is a conse¬ 
quence of the finite numbers of cold particles in the sys¬ 
tems considered in our computer simulations (see Supple¬ 
mental Material m)- However, for large systems, one 
expects that multiple clusters of cold particles should 
form throughout the system. Initially, they will grow 
due to accumulation of cold particles from the immediate 
vicinity as described above. At later stages, these clus¬ 
ters will however exhibit coarsening dynamics according 
to one or both of the following two principal mechanisms. 
Either there will be growth of larger clusters at the ex¬ 
pense of smaller clusters (Ostwald ripening) |^, or clus¬ 
ters will meet by diffusion and then coalesce. While our 
simulations do not allow us to draw any definite con¬ 
clusion with respect to the impact of Ostwald ripening, 
we can estimate the coarsening dynamics due to clus¬ 
ter coalescence. To this end, we measured the diffusion 
constant of clusters consisting of cold particles as a func¬ 
tion of their size. We hnd that the diffusion constant 
7)ci(A^coid) of saturated clusters scales in inverse propor¬ 
tion to A^coid [Fig. a- In terms of cluster radius R this im¬ 
plies that Dci{R) ^ R~^, which suggests a ‘surface diffu¬ 
sion’ mechanism for cluster diffusion [23] : Single cold par¬ 
ticles detach from the cluster surface very rapidly, then 
slowly diffuse away from the cluster surface and reattach 
to another surface site [331 [ 33 ] • This is consistent with 
our model lacking any explicit attractive interactions be¬ 
tween cold particles and their slow effective diffusion. We 
found that a cluster moving via ‘surface diffusion’ would 
grow through coalescence via R{t) ~ see Supple¬ 
mental Material m- 

Finally, we asked why-despite the lack of any explicit 
attractive interaction-there is any clustering of cold par¬ 
ticles at all. To this end we considered two cold particles 
immersed in a ‘bath’ of hot particles. As illustrated in 
Fig. ^a), in any given time interval hot particles cover 
much longer distances than cold particles, and hence will 
frequently collide with the two cold particles. In parti¬ 
cle configurations, where two cold particles are close to¬ 
gether, these collisions will tend to keep them that way, 
a phenomenon reminiscent of high-density caging of hard 
disks [35] . To quantify the degree of caging created by hot 
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FIG. 4. Diffusion constant for a saturated cluster of cold 
particles Dd in units of Dhot as a function of the total number 
of cold particles A^coid for various packing fractions ij) indicated 
in the graph; the intrinsic diffusion constant of cold particles 
was chosen as Dcoid = 0. The scaling remains unchanged 
for Dcoid > 0 (Supplemental Material [17], Fig. S7). Error 
bars indicate one standard deviation of the distribution of 
Del values obtained from ten independent simulations. The 
three straight lines indicate power law scalings Dd ~ 

K = 1.5 (dashed), k = 1.0 (solid), k = 0.5 (dotted). The best 
fit was obtained for k = 1. 


particles we determined the radial distribution function 
g{r), starting from a random configuration; see Fig.j^b). 
For Dcoid "C Dhot, gif) clearly shows that distances be¬ 
tween cold particles are likely to be shorter than when 
-Dcoid = .Dhot, where g{r) is essentially flat. We conclude 
that there is an effective attractive interaction between 
colder particles due to collisions with the surrounding 
hot particles which essentially act as a cage. This caging 
effect is fundamentally different from the well-known de¬ 
pletion interaction observed in colloidal suspensions con¬ 
taining particles with markedly different sizes IHIIMI- Ac¬ 
cording to Asakura and Oosawa [T], the effective inter¬ 
action between the larger particles is mainly an entropic 
effect: When the large particles are close together, the re¬ 
gions inaccessible to the smaller particles overlap, leading 
to an increase in the entropy of the system. 

Demixing in binary mixtures consisting of particles 
with different diffusion constants is fundamentally dif¬ 
ferent from phase separation in mixtures of Brownian 
(passive) and self-propelled (active) particles, where the 
coexisting gas and solid-like phase contain both particle 
species m- In these mixtures of active and passive par¬ 
ticles, the gas-solid phase separation occurs at very high 
Peclet numbers Pe = 3f/(2o) (where i is the persistence 
length of the active particles) and even in the absence 
of passive particles [ISlETlEH]. In contrast, in our sim¬ 
ulations the binary nature of the system is paramount 
since the differences in activity between the species drive 
the phase separation. Moreover, the binary demixing dis¬ 
cussed here occurs at Pe = 0. To compare our findings to 
the clustering at large Peclet numbers in active systems 
we also studied the case of non-zero Peclet numbers. To 
this end we performed numerical studies of binary mix¬ 
tures consisting of purely passive and active self-propelled 
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FIG. 5. (a) Illustration of caging. Two cold particles (blue, 
dark gray) are initially {t = 0) in contact and homogeneously 
surrounded by hot particles (orange, light gray). The figure 
shows two simulations at time t = 2af jDcoid with the same 
diffusion constant of cold particles Dcoid. For each time in¬ 
crement of O.Ola^/Dcoid a transparent disk with radius a is 
drawn, such that the density of disks depicts the positions 
of particles during the simulation. Top: Simulation with 
equal diffusion constants Dhot = Dcoid- Bottom: Same sim¬ 
ulation but with increased diffusion constant of hot particles 
Dhot = 10 Dcoid- (b) Radial pair distribution function g{r) for 
two cold particles surrounded by 11 hot particles at 0 = 0.13 
for three different diffusion constant ratios D = Dcoid/Dhot- 
Solid black line: g{r) of hard disks in the dilute limit ^ 0). 


particles without translational diffusion; note that trans¬ 
lational noise is not a requirement for phase separation 
in simulations of Langevin equations HHESj. For the 
model and the results see Supplemental Material HZ], 
Fig. S4 and Fig. S5. We find that demixing between cold 
and hot particles occurs even for Pe > 0. However, in¬ 
creasing the Peclet number further, demixing vanishes, 
while for very large values, we observe coexisting gas and 
solid-like phases which contain both particle species, rem¬ 
iniscent of the findings in Ref. unj. This suggests that 
demixing of passive and active particles is driven by dif¬ 
ferent mechanisms depending on the value of the Peclet 
number. 

Our predictions for binary mixtures consisting of par¬ 
ticles with different diffusion constants can be readily 
tested for mixtures of granular disks or suspensions of 
driven colloidal particles. The key ingredient for each 
system is a local mechanism that enables the particles 
to move with different diffusion constants. In the case 
of vibrated colloids, a horizontal plate is attached to a 
shaker, which induces vertical and/or horizontal oscilla¬ 
tions. By capping the plate with a lid, one can restrict 
vertical motion sufficiently to effectively confine the col¬ 
loids to two dimensions m- Systems of two disk-like 
species of identical size but with different diffusivity can 
be implemented by using particles with differently pro¬ 
filed bases [H El isa, different materials HU HU or, 
potentially, different heights. It should be ruled out that 
the collision dynamics, e.g. between particles of differ- 
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ent weight, is the cause for clustering. Driven colloidal 
suspensions represent another possible realization of our 
model [m ESI [33]. For example, consider a system con¬ 
sisting of synthetic photoactivated colloidal particles (TS] 
or non-isotropically coated Janus particles Ea iMi , to 
which a second species is added that exhibits only Brow¬ 
nian motion. Given that the persistence length of the 
active species is less than its radius, this provides two 
species with different diffusion constants. 

While our computational results show that differences 
in diffusion constant alone suffice to drive phase separa¬ 
tion in binary mixtures, the phase behavior in actual sys¬ 
tems may be even more intriguing. It will be interesting 
to explore how differences in the degree of persistence in 
the particle trajectories-especially in the regime interme¬ 
diate between the one studied here and in nni -affect the 
systems dynamics and ensuing steady states. Another 
promising route is to explore how the effective attraction 
between cold particles, mediated by caging through hot 
particles, is affected by additional interactions such as 
entropic forces mediated by differences in the size and 
shape of the particles. 
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I. IMPLEMENTATION DETAILS OF THE NUMERICAL SIMULATIONS 

In order to mimic hard particles, we use a large value for the product of spring constant and mobility kfi. We 
choose kpLO?/D\,ot = 100 with a time-discretization of df = 2 x 10“^ implying typical overlaps of \/2Dhotdt « 0.06 a 
which decay in a few time steps. We tested that the presented results are robust against reasonable changes in kpL. 
As described in the main text and in Fig. SI, the effective diffusion constant of cold particles D^s can be estimated by 
DeS = 7</'hot-Dhot) with 7 « 0.28. We simulate packing fractions 4> = (/fhot + 4’coid between 0.2 and 0.7. So even at high 
packing fractions we get D^s < 77hot/10. The relevant time scale for cluster growth is in the order of {Ll2)‘^ jDeff, 
whereas the numerical integration time step, dt, is limited by a^/Dhot- The large difference in these two time scales 
makes simulations of larger systems computationally challenging. 



0hot 


FIG. S 1. Effective diffusion constant Deft in units of Dhot of a single cold particle in a bath of hot particles at packing 
fraction <j!>hot. The diffusion constant rises roughly linear with the packing fraction (dark blue circles). To extract this density 
dependence we fit a linear function Deft = 7<('hot to the data. Note that Deft = 0 for (ft^ot ~ 0. The dashed black line shows a 
linear fit with slope 7 = 0.28. Open circles depict the diffusion constant ratio D = Dcow/Dhot at the transition points in Fig. 2, 
indicating that diffusion of cold particles is dominated by collisions with hot particles in the regime where phase separation 
occurs. For each packing fraction, the transition points are defined as the highest value of D for which Moo/Ncom ^ 0.5. 
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FIG. S 2. Finding the parameters Qatt, Ocoii, otmt through fitting. The time traces show the number M of cold particles in the 
largest cluster of the system. For each parameter set {D, (j)) a thick solid line shows the solution to Eq. (1) and thin solid lines 
depict two corresponding simulation results, (a) The parameters Oatt, Ocoii for the attachment and detachment of cold particles 
through collisions with hot particles are determined by fitting to cluster growth curves with Dcoid ~ 0. To avoid multi-site 
growth, 91 out of Aicoid = 271 cold particles are initially hexagonally close-packed. All remaining particles are distributed 
randomly. We fitted Oatt, Qcoii to time courses at a single packing fraction (j) = 0.2 and obtained ttatt = 82.5, Ocoii = 0.03. (b) 
To obtain the parameter for particle detachment due to intrinsic diffusion, aint, we use the parameters Oatt and acoii as before 
and simulate a system with D = 0.001, again at packing fraction (j) = 0.2. Choosing aint = 7.5 • 10“^ gives the correct cluster 
size at late times, (c) Similar simulations as in (a) but for different packing fractions (p. We use the parameters Oatt, Ocoii, crint 
as in (a) and (b). The solution of Eq. (1) is still in good agreement with the simulation results for different packing fractions, 
(d) Simulation as in (b) but with packing fraction cp = 0.4 and varying diffusion constant ratios D. Again we keep the fitting 
parameters fixed to the values obtained in (a) and (b). The same set of parameters describes the cluster growth at different 
densities and different D. Please note that we keep the values for all three fitting parameters aatt, cTcoII, «int fixed in every plot 
in the Supplemental Material and the main text. 
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II. RIPENING TIME EVOLUTION THROUGH COALESCENCE 

In general two ripening mechanisms exist: Coalescence of clusters and the growth of large clusters at the expense 
of smaller clusters, referred to as Ostwald-ripening [1]. Here, we ignore Ostwald-ripening in our argument and restrict 
ourselves to giving a minimal proof for ripening by only considering coalescence of clusters. Say A denotes the 
inter-particle cluster distance and 6t the time to collide with another cluster, 6t ^ Dci{R) denotes the 

diffusion constant of the cluster, which depends on the cluster radius R, Dci{R) « CojR^ (Fig. 4); Co only marginally 
depends on packing fraction thus we treat it as constant. In order to obtain a scaling for the increase of the cluster 
radius with time, we consider an arrangement of boxes of side length L; see Fig. S3 for an illustration. Each box 
contains a single cluster of radius R that covers an area fraction = ttR^/L^; similar to the one depicted in 
Fig. 1(b). This implies a mean inter-cluster distance of A « L = R^-k A collision between clusters is assumed to 
conserve the overall cluster area fraction (f>A neglecting fragmentation events during the collision process. Therefore, 
the incremental radius change during a collision is given by S In R{t) = ln(2)/2 since R + R ^ y/2R. This leads to 
-^\nR{t) = Co ln(2)/(27r). Integration yields the scaling R{t) ^ which suggests that the system ripens 

to a single cluster due to coalescence. 


L 


w 
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FIG. S 3. Illustration of the quantities used to derive the scaling argument for ripening through coalescence. The figure shows 
two compartments of an infinite system. In each compartment a cluster of radius R has formed. The cluster distance A is 
approximated by the compartment size A « L. 


HI. MODEL FOR ACTIVE AND PASSIVE PARTICLES 

Under which circumstances do binary mixtures of self-propelled particles and purely passive particles demix due to 
the same mechanisms that we observe in mixtures of particles with different diffusivities? To address this question 
we modeled mixtures of active and passive particles. Active particles perform self-propelled motion whose direction 
is subject to diffusion. Passive particles only move through collisions with hot particles. We do not add translational 
diffusion to the model. Physically, this means that we assume translational and rotational diffusion to be independent 
and translational diffusion to be negligibly small, so that the limit of small persistence lengths effectively results in 
the model described in the main text with D — 0. This assumption is only valid for systems where translational and 
rotational diffusion are not coupled. Note that in this respect the following model differs from the model described 
in [2]. However, as shown in [3], and further explained in [4], translational noise is not crucial for phase separation in 
simulations of Langevin equations. The dynamics for active and passive particles is governed by the following set of 
equations: 

There are N = Ng, + Np particles in the system, where and Np are the number of active and passive particles. F^- 
denotes the same repulsive force as described in the main text. The particles differ only in their propulsion speed: 
Vi = V for active particles and Vi = 0 for passive particles, i.e. passive particles are not propelled and move only by 
collisions with active particles. To keep particle overlaps small we choose ^ka/v = 100. The rotational noise on the 
polar vector 9i is uniform and uncorrelated: = ^SijS{t — t'), where £ denotes the persistence length of the 

active particles. The smaller the value £ the more diffusive the motion of active particles. We find that for £/a < 1 
passive particles form large clusters [Fig. S4] similar to the clusters of cold particles in mixtures of hot and cold 
particles in the main text. At these small persistence lengths we do not observe a clustering of active particles, which 
would be the case for large persistence lengths (see Fig. S5 and [2, 5]). This indicates that the effect of persistent 
motion disappears before the effect of diffusive motion becomes dominant. 
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t — 2 ■ lO'^a/v t = 2 ■ lO^a/v t — 3 ■ 10^a/i> t = 6 ■ lO^a/v 





FIG. S 4. Systems of passive and self-propelled particles at small persistence lengths behave similarly to our model. Depicted 
are snapshots for simulations with N^. = Np = 547 active (orange) and passive particles (blue), respectively. The total packing 
fraction is (p = 0.2. Different rows correspond to different persistence lengths i and different columns to different moments in 
time. For each system the particles were initially randomly distributed. For Ija > 1 only small clusters emerge. For Ija < 1, 
i.e. effectively diffusive motion, the passive particles assemble in large clusters as in simulations with particles of different 
diffusivity, where the cold particles cluster. 
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ija = 0.5 



Ija = 4 



Ija = 256 



FIG. S 5. The persistence length of the active particles determines the type of clustering. Shown are snapshots at late times of 
mixtures of active (orange) and passive (blue) particles. Initially the particles are distributed at random. For short persistence 
lengths £ < a, i.e. effectively diffusive motion, passive particles form a solid cluster. For intermediate persistence lengths i ~ a 
no large clusters emerge. For large persistence lengths £ ^ a the active particles form a jammed cluster that also contains 
passive particles. This scenario corresponds to large Peclet numbers in [2]. 



FIG. S 6. Clustering at high packing fractions. Snapshots for packing fraction 4> = 0.6 (top) and <!> = 0.7 (bottom) at different 
moments in time. Eventually a single large cluster of cold particles emerges from an initially homogeneous distribution. 
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FIG. S 7. Diffusion constants for saturated clusters Dd as in Fig. 4 in the main text. Red pentagons show the diffusion 
constant for saturated clusters in systems with packing fraction <!> = 0.5 and diffusion constant ratio D = 0.01 of cold and hot 
particles. Even at this high value for the intrinsic diffusion constant of cold particles the scaling remains Dd ~ which 

again corresponds to surface diffusion (see main text). Note that for D > 0 the values for Dd are slightly higher than for 
D = 0. This makes sense, because intrinsic diffusion of cold particles speeds up the process of surface diffusion. 
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FIG. S 8. The equal-time spatial correlation function Ce{r) = ^ -x-|=r simulations shown in 

Fig. 1 in the main text at two different moments in time. For the reference position Xi we take the position of the cluster 
particle which is closest to the center of mass of the cluster. Ce starts to decrease only at the edge of the cluster. This indicates 
that the cluster is hexagonally arranged in its entire bulk. 
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IV. VIDEO 

The video corresponds to the simulation shown in Fig. 1 in the main text. Initially both particle species are 
randomly distributed. Many shapeless clusters of cold particles emerge quickly. At longer times, a large cluster forms 
that contains most of the cold particles. This cluster remains stable. Time is measured in units of a^/Hhot- 
Parameters: N^ot = .^coid = 500, D = 10“^, (j) = 0.2. 
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